Subdiffusion and weak ergodicity breaking in the presence of a reactive boundary 
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We derive the boundary condition for a subdiffusive particle interacting with a reactive boundary 
with finite reaction rate. Molecular crowding conditions, that are found to cause subdiffusion of 
larger molecules in biological cells, are shown to effect long-tailed distributions with identical expo- 
nent for both the unbinding times from the boundary to the bulk and the rebinding times from the 
bulk. This causes a weak ergodicity breaking: typically, an individual particle either stays bound 
or remains in the bulk for very long times. We discuss why this may be beneficial for in vivo gene 
regulation by DNA-binding proteins, whose typical concentrations are nanomolar. 
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The interaction of a diffusive particle with a reactive 
boundary is of fundamental importance in interface sci- 
ence and technology, e.g., to transport in porous media 
l], interactions of proteins with artificial surfaces and 
membranes 0, or applications such as foam relaxation 
and surfactants [1 . For a Brownian particle this has been 
studied extensively, expecially concerning the question 
how bulk exchange influences the surface distribution of 
intermittently adsorbed particles Here, we derive the 
exchange dynamics with a reactive boundary of a sub- 
diffusing particle, whose unbinding and rebinding times 
in a molecular crowding environment are both shown to 
follow long-tailed distributions. We demonstrate a weak 
ergodicity breaking for the particle trajectory. 

This is of particular interest for the search of DNA- 
binding proteins for their specific binding site on DNA 
involving successive events of non-specific binding to the 
DNA and bulk excursions, such that the time spent in ci- 
ther of these events is important in the understanding of 
the various stochastic mechanisms involved in (bacterial) 
gene regulation d, @, 0] . While the generally applied as- 
sumption of Brownian diffusion of proteins works well for 
typical in vitro experiments under dilute conditions, in 
vivo the abundance of a multitude of biomacromolecules 
in the cellular cytoplasm have been shown to cause a state 
of molecular crowding: large molecules such as proteins, 
lipids, RNA molecules and ribosomes make up up to 40% 
of the cytoplasmic volume [H, 0] . In this superdense envi- 
ronment they hinder each other's motion, causing subdif- 
fusion 0, [TO, H, Hi , with a mean squared displacement 
(r 2 (i)} oc t a ; < a < 1 being a dynamic exponent 
By fluorescent methods, subdiffusion was verified for pro- 
teins in membranes with a = 0.7 [lO ]. for proteins in a 
molecular crowded in vitro environment with a = 0.75 



at higher densities [III], as well as in the cytoskeleton in 
vivo for messenger RNA of physical size ks lOOnm with 
a ks 0.75 [l2j and for dextran molecules ranging from 
lOkD to 2 MD with a in between 0.59 and 0.84 0- 
The occurrence of subdiffusion for particles with mass as 
low as lOkD was also confirmed by computer simulations 



[131 ] . The Lac repressor, a typical DNA-binding protein, 
has 141kD [lij], for which a corresponding a rs 0.73 was 
found [l3j ]. Thus, under molecular crowding conditions, 
a w 0.75 seems a fairly standard value for DNA-binding 
proteins and larger polynucleotides 11,12, II . The time 
scale over which this subdiffusion persists is not known 
precisely, but appears to be longer than minutes, so that 
the following considerations are expected to be relevant 
for genetic processes 

Esm mill- 
To derive the generalized reactive boundary condition, 
we pursue a continuous time random walk approach sim- 
ilar to Ref. [la ]: A subdiffusing particle jumps from one 
point to the next after a waiting time distributed accord- 
ing to the long-tailed probability density ip(t) ~ x" /t 1+a 
(0 < a < 1) [l7j. We start our derivation with the one- 
dimensional lattice, on which Ai is the probability to find 
the particle at lattice point i = 1, 2, 3, The probabil- 
ity of being at the reactive site (lattice point next to the 
boundary) is Ao, the notation indicating that at site 
the particle can be exchanged with the bound state with 
rate n. The balance equations then read 

dMt)/dt = I+{t)-I-(t), (la) 
dA (t)/dt = I+{t) - / "(t) - nA {t), (lb) 

and the loss from a given lattice site due to diffusion is 



JT(t) = V(«Mi(0)+ / iP(t~t')I+(t')dt' (2a) 



V(t) = *P K (t)A o {0)+ / i> K (t-t')I+(t')dt', (2b) 
Jo 

where ip K (t) = ip(t)e~ Kt . Substituting for I + from 
Eqs. |TJ, we rephrase Eqs. ([2]) in the form Jr [t) — 
f*<f>(t - t')Ai(t')dt' . The kernel $(£) is defined by 
$(w) = uip(u)/ [1 — ifj(u)} in the Laplace domain, $(u) = 
$(t)e~ ut dt [lij ]. An analogous relation holds for 
Io(t), with the kernel Q K (u) — $(u + k). For the gain 
to site i we have, assuming that the particle jumps to 
left and right equally likely, I 2 + = IjLi/2 + / l r +1 /2 and 



2 



Iq = Iq/2 + 1^/2. Note that if the particle attempts 
to jump left from site 0, it will be returned back to the 
same site. 

For the continuum limit, we introduce a new time- 
dependent quantity A (t) by J* $(t - t')A (t')dt' = 
J & K {t — t')Ao{t')dt' , corresponding to Ao(u) — 
$>(u)Ao(u)/§ K (u). Combining above results, we find 

dMt) _ r* ^ t _ t , ) A i - 1 (t')+A i+1 (t')-2A i (t>) 
Jo 



, . >. . , tit'. 

dt J G ' 2 

(3) 

In the continuum limit A(x = ai,t) = Ai(t)/a with the 
lattice spacing a [lil ] , this equation for i > 1 yields 



dA(x, t) a 



Of 



■2 I *< ! 



0^*. (4) 



In the long time limit w <C 1, ~ 1 — (ut)", and 



$(u) 



V~~ Q to leading order. With AT Q = a 2 /[2r Q ] 



and the fractional Ricmann-Liouville operator 

1 d f* A(x,t') 



oD 



l-a 
t 



A(x,t) 



T(a)dtJ (t-f) 



l\l-a 



dt\ 



(5) 



Eq. (HJ) is equivalent to the fractional diffusion equation 



dA(x,t)/dt = K aa Dl~ a d 2 A{x,t)/dx 2 for x > [14 1. 
Similarly, Eq. (|lb[) can be recast into the form 



d . , „ 



(t-i')A(t')^' 
^(tfJ-AKf) 



(6) 



In the continuum limit, we recover the expression 



S(t)A (0) 



V(t-t')A(Q,t')dt' 



t') 



dA{x, f) 



dx 



dt' (7) 



with <3>(m) = o(u+k)$(«)/$ k (u). -4 (0) is 1, if the parti- 
cle is initially released at site 0, and otherwise. The re- 
action rate at the boundary is Jreact = aK Jn v&k 1 ^) — 
t')A(0, t')dt' , and the right hand side of Eq. represents 
the flux into a; = from positive x. We expand Eq. at 
?i = in Laplace space (note that u <C K [19fl) producing 
the sought for reactive boundary condition 



-A(0) + fco^r"A(0,i) (8) 



x=0 



for the subdiffusive particle. This is one of the main re- 
sults of this work. We defined k = 2nK a / [a$ K (u = 0)] ~ 
ok q , using that kt — > in the continuum limit [l9| . 

The Berg-von Hippel model maps the bind- 
ing/unbinding dynamics of a DNA-binding protein 
to/from the DNA surface onto a cylinder of radius R\ 



laced along the z-axis in cylindrical coordinates (r, 6, z) 
For particles subdiffusing in the space r > R\ with 
density P the boundary condition (jHJ) generalizes to 



2nR 1 K aQ Dt a dP(r,t)/dr\ r=Ri 



= -P /L + k on D^ a P\ 



r=R! ■ 



(9) 



where k on = 2irRik, and we have assumed rotational and 
translational symmetry around and along the z-axis. L 
is the length of the cylinder along the z-axis, and Pq = 1 
if at t = the particle is at the boundary, and Pq = 
otherwise. The rate of reaction with the cylinder per 
length along the z-axis is j Icac t{t) = k on0 Dl~ a P\ r=Rl , 
and the fractional diffusion equation ((4]) is replaced by 



dP x _ a l d ( BP 



(10) 



Consider now the situation when the particle is bound 
at the boundary at t = 0. In a crowded environment 
the full escape of the particle from the boundary consists 
of two steps: unbinding with rate K g, returning it to 
the exchange site 0, and then avoiding to rebind to the 
boundary such that the particle is unbound at the mo- 
ment when the environment allows for jumping to site 1. 
If the particle is bound at the time when the environment 
allows a jump (in principle), then the process needs to 
start over, etc. Thus, the waiting time distribution for 
the full escape can be written as 



(11) 



Punb(i) = R(t) + R(t- t')Tl(t')dt' 
Jo 

+ f R(t- t') f U(f - t")n{t")dt"dt' 
Jo Jo 



With R(t) = V(«)[l-Pbou n d(*)] andn(t) = ^ (Abound (*)• 

Here Pbound(^) is the probability that the particle is 
bound at time t, given that it was bound at t = and 
that the environment has not yet opened up to allow a 
jump to i = 1. This yields 



Pbound (^) 



Koff 



-(K+Koff)* 



K + K„ff K + K S 



(12) 



The waiting time ip(t) between opening events in the en- 
vironment, allowing a jump to i = 1, is given by 



(13) 



where the factor 1/2 accounts for the fact that the parti- 
cle jumps to site 1 only with a probability 1/2 (and stays 
at with probability 1/2). In Laplace space, ip(t) can 
be expressed in closed form through a geometric series, 
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the bound state p b 



FIG. 1: The distribution 5 a , Eq. ([18)) , for various /3, with 
a = 0.75. In all divergence at p = and 1 is observed. 

The points are results from a stochastic simulation [2^|. 



^(«) = [V»(«)/2]/[l-^(u)/2] 
the Laplace domain, 



1 - 2(W)". Similarly, in 



punb(w) - $(«) - n(«) + $(«) - n(u)]n(«) 
+[f( il )-nM]n(«f + .-= n i ; —A ^ , (14 ) 



where 11(f) in Laplace space assumes the exact form 

n(u) = — - — 4>{u)+ KoS i>{u + K + KoS ). (15) 

K + K s K + KoS 

Collecting the results, we obtain for p un b at small u (ljj 



Punb(u) ~ 1 



K + K S 



K oe (K + K s) a 



l-u a /k s, (16) 



where in the continuum limit (k q s -C k) we have k Q s ~ 
Kos/k 1 ^". The unbinding times then are distributed 
according to the power-law p un b — 1/ (fc ffi 1+Q! )- This 
is a central finding of this work: The crowded environ- 
ment impeding the desorption to the bulk translates the 
a priori exponential distribution of unbinding times to a 
power-law [20]. Once arrived at site i = 1, the particle 
subdiffuses in the bulk. We consider here the cylindrical 
case governed by Eq. (fl"0|) . With initial condition P = 1 
and a reflecting boundary condition &tr = R 2 jll[ , an an- 
alytic result can be obtained in terms of modified Bessel 
functions, see Ref. [22| for details. A systematic expan- 
sion for small u leads to the result [2211 



Prob(u) ~ 1 - Su a /k 



(17) 



with the cylindrical cross-section S = n(R 2 — R 2 ). The 
form p re b(t) ~ f~ 1_Q is typical for subdiffusion [14 1. 

Both unbinding to the volume and returning to the 
reactive boundary follow power-law forms with identi- 
cal asymptotic behavior ~ <~ 1-Q . The lack of a char- 
acteristic time scale separating micro- and macroscopic 
events gives rise to weak ergodicity breaking 23j. As 
shown in Ref. 24], the time-averaged probability in 



gle trajectory, 

Sot(koa/(Sk s] 
el-function 24], 



Pbo 



ound = hm t _oo <boundA f° r a sin ~ 

has the distribution V (p bound ) = 



Pbound ) ' w hh the Lamperti-generalized 



5a(P,P) 



7T" 1 sin( 7 ra)/3p Q - 1 (l - p)"" 1 
/3 2 (1 - p) 2a + p 2a + 2/3(1 - p) a p a cosTra' 

(18) 

Note that V is normalized, f Q 'P(Pb und) c ^Pbound = 1> an d 
valid in the long t limit. It is independent of t and in that 
sense an equilibrium is attained. However, while in the 
Brownian limit a = 1, ergodicity and a sharply peaked 
behavior for V are recovered, the very distinct behav- 
ior of V for a < 1 is displayed in Fig. [T] for a = 0.75: 
as function of p = p bound , the distribution peaks at 
and 1, with a smaller maximum in between. Thus, in 
a single trajectory a particle is typically either bound 
or unbound, independently of the duration of the tra- 
jectory. This nonergodic behavior is imposed on the 
system by the probability ip(t')dt' ~ t~ a of never 
moving, that decays very slowly. The smaller the cross 
section S, the more likely is it to find the particle in a 
bound state, as it should be. The behavior of V therefore 
contrasts the ensemble average over many trajectories, 
(Pbound) = (1 + Skos/kon) \ corresponding to the form 

Abound ) = HPbou 



- k on /[k on + Sk s]) [24j- This can 
be understood as follows. For an ensemble of particles, 
kon/koS defines the nonspecific binding constant K I1Sl 
equal to the ratio N oouna \/ (S N unoouna \) of bound and un- 
bound particles normalized by the cross section [?| ■ Then 
1/(1 + Skos/kon) is the ensemble probability that a par- 
ticle is bound. Weak ergodicity breaking is thus relevant 
for systems with few particles of a given species. 

Transcription factors (TFs), DNA-binding proteins 
regulating the transcription of a specific gene, occur at 
very small numbers (a few to some hundred per cell 26]), 
and in many cases it is essential for the stability of ge- 
netic circuits that a TF is always bound at some oper- 
ator site on the DNA HI]. While the random mo- 
tion of the TFs in most in vitro experiments is Brow- 
nian, molecular crowding in vivo causes subdiffusion of 
TFs. This would have interesting consequences for gene 
regulation. Namely, due to the weakly ergodic behavior 
demonstrated here, TFs will typically stay close to their 
binding site with a diverging characteristic time scale, 
such that unbinding and escape to the volume is greatly 
reduced. The price to pay is that once a TF escapes to 
the bulk, its return is also affected by an infinite average 
time. Moreover, there exists a large class of TFs, such 
as the well-studied Lac and bacteriophage A repressors 
in E. coli (28|, whose specific binding site is located im- 
mediately adjacent to their coding region. Biochemical 
production occurs likely within a few tens of nm from 
the coding region [2§], and therefore from the targeted 
binding site. The weak ergodicity breaking thus keeps 
those TFs within a small volume around their complete 
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biochemical cycle, very likely leading to a significant in- 
crease in the stability of the regulation of that partic- 
ular gene. Subdiffusion caused by molecular crowding 
could therefore be very beneficial for living cells, allow- 
ing them to maintain the concentrations of even vital 
TFs at nanomolar levels. This may significantly impact 
our current picture of gene regulation in vivo and pose 
the need to perform experiments much closer to the cel- 
lular crowding conditions in order to obtain meaningful 
information for the in vivo situation. 

We derived the generalized reactive boundary condi- 
tion for the interaction of a subdiffusive particle with 
a boundary and showed that in the molecular crowd- 
ing scenario the distribution of unbinding times becomes 
long-tailed, with the same exponent as the distribution 
of return times to the boundary. This gives rise to weak 
crgodicity breaking, relevant for systems with small num- 
bers of diffusing particles. Apart from gene regulation, 
these effects will impact cellular processes in more gen- 
eral, such as the interactions of biopolymers with mem- 
brane proteins, or the exchange of shorter DNA and RNA 
chains across cellular membranes. Moreover, they will af- 
fect trapping phenomena in the vicinity of soft interfaces 
in more general, e.g., the exchange dynamics from ion 
clouds in the vicinity of charged or polarized membranes. 
It should be very interesting to explore these effects by 
single particle tracking under molecular crowding condi- 
tions using fluorescent labelling techniques. 
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